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Abstract. Let B H be a fractional Brownian motion with Hurst index < H < 1/2. 

■ In this paper we study the generalized quadratic covariation [f(B ),B ]^ ' defined by 

-I— > . 



oo 



where the limit is uniform in probability and x h-> f(x) is a deterministic function. We 
£Nj , construct a Banach space Jff of measurable functions such that the generalized quadratic 

^ ' covariation exists in L 2 and the Bouleau-Yor identity takes the form 

O . [f(B H ), B B ] W = _ / f( x )J? H (dx, t) 

provided / £ where J? H (x,t) is the weighted local time of B H . This allows us 
to write the fractional Ito formula for absolutely continuous functions with derivative 
belonging to Jtf. These are also extended to the time-dependent case. 

1. Introduction 

Given H £ (0, 1), a fractional Brownian motion (fBm) with Hurst index H is a mean 
zero Gaussian process B H = {Bf, 0<t<T} such that 

E [B?B?] = i [t 2H + s 2H -\t- s\ 2H ] 

for all t, s ^ 0. For = 1/2, coincides with the standard Brownian motion B. B H is 
neither a semimartingale nor a Markov process unless H = 1/2, so many of the powerful 
techniques from stochastic analysis are not available when dealing with B H . As a Gaussian 
process, one can construct the stochastic calculus of variations with respect to B H . Some 
surveys and complete literatures for fBm could be found in Biagini et al [2] , Decreusefond 
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and Ustiinel [BJ, Gradinaru et al [3ESJ, Hu [IS], Mishura [H] and Nualart [23]. It is well- 
known that the usual quadratic variation [B H , B H ]t = for 2H > 1 and [B H , B H ]t = oo 
for 2H < 1, where 

[B H ,B H ] t = lim± f(Bf +£ -Bffds 

eiU £ Jq 

in probability. Clearly, we have also 

3=1 

where the limit is uniform in probability. This is inconvenience to some studies and 
applications for fBm. We need to find a substitution tool. Recently, Gradinaru et al [2] 
(see also |15] and the references therein) have introduced some substitution tools and 
studied some fine problems. They introduced firstly an ltd formula with respect to a 
symmetric-Stratonovich integral, which is closer to the spirit of Riemann sums limits, and 
defined a class of high order integrals having an interest by themselves. On the other hand, 
inspired by Gradinaru-Nourdin [T2J [13] and Nourdin et al [211 122] j as the substitution tool 
of the quadratic variation, Yan et al |27j considered the generalized quadratic covariation, 
and proved its existence for ^ < H < 1 (Thanks to the suggestions of some Scholars we 
use the present appellation). 

Definition 1.1. Let < H < 1 and let f be a measurable function on R. The limit 

(1.1) Hm§ f {f(B* +£ ) - f(B*)} {B? +£ - Bf) S ™- l d S 

is called the generalized quadratic covariation of f{B H ) and B H , denoted by [f(B H ), B H ][ W 
provided the limit exists uniformly in probability. 

In particular, we have 

[B H ,B H ]l W) =t 2H 

for all < H < 1. If H = i, the generalized quadratic covariation coincides with the 
usual quadratic covariation of Brownian motion B. For \ < H < 1, Yan et al [28] showed 
the generalized quadratic covariation can also be defined as 

(1.2) [f(B*),B*\M=2H hm V (A;) 2 *" 1 {/«)-/«,)}«- <J, 

|7r„| — »0 

tj fc7Tn 

provided the limit exists uniformly in probability, where ir n = {0 = to < t\ < • • • < t n = t} 
denotes an arbitrary partition of the interval [0, t] with \\Tr n \\ = supjitj — tj-i) — > 0, and 
Aj = t .jj.. , j = 1,2,... ,n. Moreover, by applying the time reversal Bp = B^_ t on 
[0, T] and the integral 

/ f(x)^ H (dx,t), 

Yan et al [28] constructed a Banach space Mh of measurable functions such that the 
generalized quadratic covariation [f(B H ), B H ][ W ^ exists in I? if / G B#, where 

Sf H (x, t) = 2H f 5{B* - x)s 2H - 1 ds 
Jo 
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is the weighted local time of fBm B H . However, when < H < | the method used in Yan 
et al [271 EH] is inefficacy. In the present paper, we shall consider the generalized quadratic 
covariation with < H < i. Our start point is to consider the decomposition 

ft 



(1.3) 



£ Jo 



i / f{Bf +£ )(B^ £ -B?)ds 2H --^ ! f(B»)(B» +£ -B*)ds™ 
£ Jo £ Jo 



Clearly, if the modulus in expression (|1.3p is ~, the decomposition is meaningless in general. 
For example, for f(x) = x we have 

\ f E [B?(B? +£ - Bf)] ds™ = \f\[(B + e? H - - e™] ds™ 
£ Jo £ Jo 1 



-co, 



as e I 0. However, 

~2H I \ EB ^ ( B ^+e - B ^ )\ ds2H 

£ Jo 

as e I 0. Thus, for < H < \ we can consider the decomposition (|1.3|) . By estimating the 
two terms of the right hand side in the decomposition (jl.3p . respectively, we can construct 
a Banach space of measurable functions / on R such that ||/||jsf < oo, where 





V^(T-s) 1 -"' 

We show that generalized quadratic covariation [f(B H ), B H ][ W ^ exists in I? for all t G 
[0, T] if / G J^. This allows us to write Ito's formula for absolutely continuous functions 
with derivative belonging to Jff and to give the Bouleau-Yor identity. It is important to 
note that the decomposition (jl.3j) is inefficacy for | < H < 1. 

This paper is organized as follows. In Section [2] we present some preliminaries for fBm. 
In Section [31 we establish some technical estimates associated with fractional Brownian 
motion with < H < | • In Section 01 we prove the existence of the generalized quadratic 
covariation. We construct the Banach space Jf? such that the generalized quadratic co- 
variation [f(B H ), B H ]( W ^ exists in L 2 for / 6 Jtif. As an application we show that the ltd 
type formula (Follmer-Protter-Shiryayev's formula) 

F(B H ) = F(0) + jf f(B?)dB? + \ [f(B H ), B H ] ™ 

holds, where F is an absolutely continuous function with the derivative F' = f G 3if. In 
Section [5j we introduce the integral of the form 

(1.4) / f(x)J? H (dx,t), 
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where x i— >■ f(x) is a deterministic function. We show that the integral (|1.4[) exists in L 2 , 
and the Bouleau-Yor identity takes the form 



\f{B B ),B H ]W = - I f(x)J? H (dx,t) 



provided / G Jff. Moreover, by applying the integral (|1.4|) we show that (jl.ip and (jl.2|) 
coincide for < H < 5 when / G In Section [6l we consider the time-dependent case, 
and define the local time of B H with < H < \ on a continuous curve. 

2. Preliminaries 

In this section, we briefly recall some basic definitions and results of fBm. For more 
aspects on these material we refer to Biagini et al [2], Hu |18] . Mishura |19j . Nualart [23] 
and the references therein. Throughout this paper we assume that < H < \ is arbitrary 
but fixed and let B H = {B? ,0 < t < T} be a one-dimensional fBm with Hurst index 
H defined on (O, J 7 , P). Let (S)* be the Hida space of stochastic distributions and let o 
denote the Wick product on (5)*. Then 1 1-> B^ is differentiable in (5)*. Denote 

We call W(^) the fractional white noise. For u : R + — > (5)*, in a white noise setting we 
define its Wick-Ito-Skorohod (WIS) stochastic integral with respect to B H by 

(2.1) I u s dB? := f u s oW^ H) ds, 

Jo Jo 

whenever the last integral exists as an integral in (5)*. We call these fractional ltd in- 
tegrals, because these integrals share some properties of the classical ltd integral. The 
integral is closed in L 2 , and moreover, for any / G C 2,1 (]R x [0, +00)) the follwing ltd type 
formula holds: 

f(B t H , t) = /(0, 0)+ j* ^f(Bf, s)dB» 

(2.2) + j* §- g f(B?, s)ds + H jT ^f(B?, s)s 2H - 1 ds. 

The fBm B H has a local time C H (x, t) continuous in (x, t) G R x [0, 00) which satisfies the 
occupation formula (see Geman-Horowitz |llj ) 

(2.3) f 4>(B*,s)ds= [ dx [ <f>(x, s)C H (x,ds) 

Jo Jr Jo 

for every continuous and bounded function cj>(x, t) : R x R + — > R, and such that 

rt 1 
£ H (x,t)= / 6(Bf-x)ds = ]im—X(s£[0,t],\B^-x\<e), 
Jo 40 2e 

where A denotes Lebesgue measure and 5(x) is the Dirac delta function. Define the so-call 
weighted local time J£ H (x,t) of B H at x as follows 



Jf H (x,t) 



2H f s 2H - 1 C H (x,ds) = 2H f 5{B* -x)s 2H - l ds. 
Jo Jo 
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Then the Tanaka formula 
(2.4) 



B? - x\ = \x\ + [ sign(£f - x)dBf + £? H (x, t) 
Jo 



holds. 

For H £ (0, 1) we define the operator M on L 2 (R) as follows (see Chapter 4 in Biagini 
et al [2J and Elliott- Van der Hoek [8]): 

/3h d f (s — x) 



Mf(x) 



H -\dx 



\S — X 2 



f(s)ds, f G L'(R), 



where (3h is a normalizing constant. In particular, for H = i we have Mf(x) = /(x), and 



for < H < | we have 



M/(x) = Pe I f ^- S )-J ix) ds. 



S 2 



-11 



As an example let us recall Ml[ a () ](x), i.e., M/ when / is the indicator function of an 
interval [a,b] with a < b. By Elliott- Van der Hoek [8], Ml[„^(i) can be calculated 
explicitly as 



(2.5) Ml [aM (x) 



b — x 



yjT{2H + 1) sin(vr^) 
2r( J ff + i)cos(f(i? + i)) I |6-x| 



a — x 



O — X 2 



-iJ 



By using the operator M we can give the relation between fractional and classical white 
noise (see Chapter 4 in Biagini et al [2]) 



which leads to 



u t dBf 



M (ul [0jT] ) t SB t , 



where u is an adapted process and J R vt5B t denotes the Skorohod integral with respect to 
Brownian motion B defined by 



VtSBt 



v t o W t dt. 



( H) 

Let D t denotes the Hida-Malliavin derivative with respect to B n . In the classical case 
{H = 1/2) we use the notation Dt for the corresponding Hida-Malliavin derivative (for 
further details, see Nualart [23] and Biagini et al [2]). We have 



D t F = MD { t H) F 



and 

(2.6) E 
for F £ L 2 (P). 



F / u s dBf 
Jo 



E 



{Mul [m ) s {MD^F)ds 
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3. Some basic estimates 
In this section we will establish some technical estimates associated with fractional 
Brownian motion with < H < |. For simplicity throughout this paper we let C stand 
for a positive constant depending only on the subscripts and its value may be different in 
different appearance, and this assumption is also adaptable to c. 

Lemma 3.1. For all t, s S [0, T], t > s and < H < 1 we have 

(3.1) 1(2 - 2 H )s 2H (t - s) 2H < t 2H s 2H -y 2 < 2s 2H (t - s) 2H , 

where y. = E(B^B^). 

By the local nondeterminacy of fBm we can prove the lemma. Here, we shall use an 
elementary method to prove it. We shall use the following inequalities: 

(3.2) (1 + x) a < 1 + (2 a - l)x a 

(3.3) (2 - 2 a )x a {\ - x) a < (1 - x) a - (1 - x a ) < x a (l - x) a 

with < x,a < 1. The inequality (|3.2p is a calculus exercise, and it is stronger than the 
well known (Bernoulli) inequality 

(l + x) a < l + ax a < l + x a , 

because 2 a — 1 < a for all < a < 1. The inequalities (|3.3p are the improvement of the 
classical inequality 

1 ~x a < (l-x) a . 
The right inequality in (j3.3[) follows from the fact 

(l-x) a (l-x a ) <l-x a . 

For the left inequality in (I3.3p . by (I3.2p we have 

l = (l-i + i) a <(l-x)"V/ + (2 a - 1) [(1 - x) a A x a ] 

for < x < 1, where x V y = max{x, y} and x A y = minjx, y}, which deduces 

(1 - x) a - (1 - x a ) > (2 - 2 Q )(1 - x) a A x a 
> {2-2 a ){l-x) a x a . 

Proof of (|3,ip . Take s = xt, < x < 1. Then we can rewrite p rjS := t 2H s 2H — fi 2 as 

Pr , s = t± H L 2H -\[l + x 2H -{l-x) 2H ] 2 \ 
= t 4H G(x). 
In order to show the lemma we claim that 

(3.4) 1(2 - 2 H )x 2H (l - x) 2H < G{x) < 2x 2H {\ - x) 2H 
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for all x G [0, 1]. We have 

G{x) = x 2H -\[l + x 2H -{l-x) 2H f 

= I {2x H - (1 + x 2H - (1 - x) 2H ) } {2x H + (1 + x 2H - (1 - z) 2 ") } 



{(1 _ x )2^ _ (1 _ x ")2} | 2x H + X 2H + 1 _ (1 _ x) 2H| 

{(1 - x) H - (1 - x H )} {{l-x) H + l- x H ) {2x H + X 2H + 1-(1- x) 2H } . 



4 

Thus, (13.41) follows from (13.31) and the facts 



(1 - x) H < (1 - a;) H + (1 - x H ) < 2(1 - x) H , 
2x H < 2x H + x 2H + 1-(1 - x) 2H < 4x H . 

This completes the proof. □ 

Lemma 3.2. For all t, s € [0, T], t > s and < iJ < | we /iave 

(3.5) i(t- S )^<t^-/i<(t- S r, 



and 



(3.6) -(2 - 2")(-)- (t - ,)- -(-)-!/ - , s 

w/iere /i = E(B^Bf). 
Proof. The inequalities (13. 5p follow from 



t 2H - p = t 2H - L - {t 2H + s 2H - (t - s) 2H ) 

1 /.9W ou\ 1 



L{f»-s 2H ) + t{t-s? H . 



In order to show that (13.61). we have 



s ™ _ ^ = S 2H _ 1 + g 2H _ (t _ s) 2tf) 

W(i-*r-(i-(fr 



Thus, the inequalities f|3.6|) follow from (|3.3p . This completes the proof. □ 
Lemma 3.3. For < H < ^ we /wroe 

(3.7) \E [{B? - B?)(B» - B?)] \ < cj* ~ ^ " ^ 



2H 



for all < s' < t' < s < t. 
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Moreover, the estimate (|3.7p holds also for all < s' < s < t' < t. In fact we have 





(t' 


-s) 4H 


- (t' - 






?H(1/ - S ') 2H 
' V 15 ) i 


(4 


t ) [t 


- s ) 


< ft 
S (t - 


S) {t - 


J\2H 

- S ) ! 




(s~ 


s') 2H (t' 


-s) 2H 




-s') 2H {t- 


-s) 2H , 






(*- 


- s') 2H 


= {(*" 


-s) + (s 


-s')} 2H <(t- 










_ it- 


s) 2H (t' - 


- s) 2H + {s- s r 


)2H( t ,_ s) 2H 












(f - s) 2H 










it- 


- s) 2H (t' 


- s') 2H 












(t> - s) 


2H ' 





which gives 

\E [{Bf - Bf ){Bf - Bf)] | = i {|i - + |* - f I 2 * - |t - f I 2 * - |* - 

< o (*-«) 2g (* , -^) 2g 

(tf - s) 2H 

Proof of (|3.7p . For < s' < t' < s < i < T we define the function x i— >■ G s j(x) on [s',t'] 
by 

G s>t (x) = (s - x) 2H - (t - x) 2H . 
Thanks to mean value theorem, we see that there are £ G ( s '>t') and r\ E such that 

2E [(Bf - Bf)(Bf - Bf)] = G., t (0 - G., t (*0 

= 2ff (f - s') [(* - 2H ~ l ~ (s - O 2 "' 1 } 



2H{2H - l){t' - s')(t - s) (t? - i) 2H - 2 < 0, 



which gives 



(3-8) \E [{Bf - Bf){Bf - Bf,)] | < ■ 

On the other hand, noting that 

(i-s)"(t'-s') H 



we see that 



\E[(Bf - Bf)(Bf - Bf)]\ ( \E [{Bf - Bf)(Bf - Bf)] \ 



{t- s ) H {t> - S ') H -y (t - s) H {f - S ') H 

for all a G [0, 1]. Combining this with (I3.8j) . we get 
\E[{Bf - Bf)(Bf - Bf)}\ < 



H tz>H\ ( ~dH oH ,^(t-s)^ H+a (t'-sT- a)H+C 



and the lemma follows by taking a = H/(l — H). □ 
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Lemma 3.4. For < H < \ we have 

\E[B«(B t H -B?)]\<(t-s) 2H , 
\E[B?(B?-B?)]\<(s-rf H , 
\E[B?(B»-B?)]\<(t-s) 2H 

for all t > s > r > 0. 

Let tp(x,y) be the density function of (B^,B^) (s > r > 0). That is 

/ \ 1 I 1 / 9 9 9W 9\ I 

where /i = E(B^B^) and p 2 = r 2H s 2H - i? . 

Lemma 3.5. Let f G C 1 (M) admit compact support. Then we have 

\E [f{B*)f{B*)] | < rg( ^ )2g (J5 [\f(B? )| 2 ] £ )| 2 ]) 1/a 
/or all s > r > and Q < H <\. 
Proof. Elementary calculation shows that 

f 2 {y){x - ^j W y) 2 (p{x,y)dxdy 

which implies that 
1 



P 4 



\f{x)f(y){s 2H y-^x)(r 2H x - py)\ip(x, y)dxdy 



< 



r S (T? \\t(T3H\\2-] r* ntr-DH\m\V i 



P 2 



E[\f{B*)\*\E[\f{B*)\*\y 



± ^11^ & 0/(^)| 2 ] E [|/(^)| 2 ]) 1/2 



by Lemma 13. II It follows that 

\E[f'(B?)f(B?)]\ = \ [ f{x)f(y) 1 f-^x,y)dxdy\ 

J R 2 dxdy 

f( x )f(y) (s 2H V - iix)[r 2H x - ny) + j ip{x, y)dxdy\ 

^ r nl H l H rf H {E[\f{B?)\ 2 ]E[\f { B?)\ 2 ]) l '\ 
This completes the proof. 

Lemma 3.6. Let f £ C 2 (M) admit compact support. Then we have 

\E [f"{B»)f{B?)\ | < (E [\f{Bf )| 2 ] E [|/(^)| 2 ]) 1/2 



for all s > r > and < H < \ 
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Proof. A straightforward calculation shows that 



f 2 (y)( x ~ -§Hy) 4 ^(x,y)dxdy = ~4H f(y) r— e'^dy, 



which deduces 



-4 / f(x)f(y){r 2H x-fiy) 2 (p(x,y)dxdy 



P 4 



C H 

(s — r) 5 



< T, " jE[\f(B»)\*}E[\f(B*)\i] 



by Cauchy's inequality and Lemma 13. 11 It follows that 

\E\f»(B?)f(B?)]\ = \J f(x)f(y)-jj^<p(x,y)dxdy\ 



< 



f 1 r 2H ) 
f(x)f(y) <i-^( r2Hx ~ VV) 2 ~ — J f(x,y)dxdy\ 

C H ,1/2 



E[\f(B^]E[\f{B^]Y'\ 



(s - r) 2H 

This completes the proof. □ 
4. Existence of the generalized quadratic covariation 

In this section, for < H < \ we study the existence of the generalized quadratic 
covariation. Denote 

J £ (f,t) := f - f(B? )} {Bf +£ - Bf )ds 2H 

£ Jo 

for e > and t > 0. Recall that the generalized quadratic covariation [f(B H ), B H ]^ W ' is 
defined as 

(4.1) \f(B H ),B H ]W :=limJ £ (/,t), 

e4-0 

provided the limit exists uniformly in probability. Clearly, we have (see, for example, 
Klein and Gine 1 16 ) 



(4.2) [B H ,B H ]W=t 2H 
for all t > 0. In fact, one can easily prove that 

2 



E 



[\ B « +£ -B?) 2 ds-t 
e Jo 

= i f f E " ^) 2 «, " Bff] ds 2H dr 2H - t* H 

e Jo Jo 



for t > 0, as e ! 0. 
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Consider the decomposition 
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£ Jo 

(43) = -i f /«j(^f +£ - ^ ^ - i r /(^f)« £ - 

£ Jo £ Jo 

= I+(f,t)-I-(f,t), 

and define the set Jt? = {/ : measurable functions on R such that < oo}, where 



dxds 



V2^{T- S y- H ' 

Then, Jf 7 is a Banach space and the set <? of elementary functions of the form 

i 

is dense in J4f, where {xi,0 < i < 1} is an finite sequence of real numbers such that 
Xi < Xi + \. Moreover, Jtf? contains the sets J^, 7 > 2, of measurable functions / such that 



Our main object of this section is to explain and prove the following theorem. 

Theorem 4.1. Let < H < \ and f G Jtf. Then the generalized quadratic covariation 
[f(B H ),B H ]W exists and 



(4.4) E [f(B H ),B H ]^ 



We split the proof into several lemmas, and for simplicity throughout this paper we let 
T= 1. 

Lemma 4.1. Let < H < \ and let f be an infinitely differentiable function with compact 
support. We then have 

(4-5) E\l-(f,t)\ 2 <C H \\f\\^, 
(4-6) E\lt(f,t)\ 2 <C H \\f\\^ 



for all < e < 1. 
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Proof. We need only to obtain the first estimate. It follows from (j2.6f) that 



E [f(Bf)f(Bf)(Bf +£ - Bf )« e - Bf)] 

rr+e 

= E f(Bf)f(Bf)(Bf +£ - Bf) / dBj 

Jr 

= E [ Ml [r . r+£] {l)MD\ H) f(Bf)f(Bf)(Bf +£ - Bf)dl 
Jr 

Ml [r>r+e] (l)Ml [0;S] (l)E [f'(Bf)f(Bf)(Bf +£ - Bf)] dl 

Ml [r , r+E] (l)Ml [Q , r] (l)E [f(Bf)f(Bf)(Bf +£ - Bf)] dl 



+ / Ml [ryr+£] (l)Ml [SyS+£] (l)E[f(Bf)f(Bf)]dl 
Jr 

= E [Bf(Bf +£ - Bf)] E [f'(Bf)f(Bf)(Bf +£ - Bf)] 

+ E [Bf(Bf +£ - Bf)] E [f(Bf)f'(Bf)(Bf +£ - Bf)] 
+ E [(Bf +£ - Bf)(Bf +£ - Bf)] E [f(Bf)f(Bf)] 
= V £ (s, r, 1) + * e (a, r, 2) + * e (s, r, 3). 



In order to end the proof we claim now that 



(4.7) 



T 4H 



f f y £ (s,r,k)ds 2H dr 2B 
Jo Jo 



<C H \\f\\ 



2 



k — 1 ; 2, . 3 , 



for all e > small enough. Some elementary calculus can show that, for all < e < 1 



r E[\ f {Bft]s™-Hs r e - 

Je JO r 



dr 



2H (s — £ — r) 2H 
fE[\f{Bf)?] s^ds H- 

Je JO r 

s 2H - l E[\f(Bf)\ 2 ] ds 
dr 



dr 



~ 2H (s — e — r) 2H 
1 dr 



o x^il-x) 211 



dx , 



E[\f(Bf)\ 2 ]s 2 ^ds , rl _ 2H[r + £ _ s)2H 



< J X E[\f(Bf)\ 2 ] da J' 



E[\f(Bf)\ 2 ] ds 



dr 



,2H 



< J^E[\f{Bf)\ 2 ] cfaQf 



r 2 iH (r + e — s 

: dx 

x 2 - m (x - 1) 2H 
+°° dx 
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and 



E[\f(Bf)f]s 2H -'ds 



r dr 



o {r + e-s) 2H 
E[\f(B?)f] s^ds 



x 2H ' l dx 

Y2H 



< / E{\f(B?)\ 2 )ds 



o 0- + x) 
:;,/ -' < r°° x H ~ l dx 



(4.8) 



io (e-s) H \J (1 + x) 

Jo Jr V2tt(1 - s) H 



2H 



where the estimate (14. 8p follows from the monotonicity of the function 



e i-)- 



e g 2H-l 



(e-s) 



e 2s' 2H ds 



with e G [0, 1]. It follows that 



l pi 




2H,2H 



■$f e (s,r, 3)ds 2H dr 



< 



10 JO 

H rl 



J J \E [(Bf +£ - B*)(B& - Bf)] | 

• {E [/ 2 (Bf )] + E [/ 2 (Bf )] } {srf^dsdr 

4h t t \ E [«, " B r)(B? +£ ~ Bf)] \E [/ 2 (Bf )] (sr) 
£ Jo Jo 



2H-1 



dsdr 



<H I E[\f{B»)\*]s™- l ds 



dr 



A-2H 



(s — e — r) 



2H 



+ H^E[\f(Bf)\ 3 ] s 2H ^ds 



dr 



H\\2~\ 2H-1 



+ H / E[\f(B;')r\ , 
Jo 



ds 



,l-2Hl 



s r 2H-l dr 



\2H 



\2H 



io yr + e-8) 
<C H \\f\\%> 

for all < e < 1. 

Now, let us obtain the estimate (|4.7|) for k = 1. By (j2.6[) we see that 
r, 1) = B [Bf (Bf +£ - Bf )] B [/'(Bf )/(Bf )(Bf h£ - Bf )] 

= B [Bf (flg. - Bf )] B [flf (Aft. - Bf )] B [r (Bf )] 

+ ^ [flf (5^ - Bf )] B - Bf )] B [f'(Bf)f'(Bf) 

= tt e (s,r,l,l) + tt e (a,r,l,2). 
Together Lemma 13.51 Lemma 13.61 Lemma 13.41 and the fact 



(4.9) 



£ [/W)] = / Z 2 (*) 



2irr H 



-e zr 2H dx 



H 
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with s > r > lead to 



1 

TUT 



f i\ £ {s,r,l,l)ds 2H dr 2H < f f\Elf"{B»)f{B?)]\ds 2H dr w 
Jo Jo Jo Jo 



<C H 



t r s 



JO 

t 



(s - r) 2H 



< C H I E[f 2 (B?)]ds 2H 



E\f(B?)f(B?)\ds 2H dr 

s H/2 



o (s — r ) 2H r H / 2 



dr 



2H 



<c H 



and 



e 



1 

AH 



f fMs,r,l,2)ds™dr™ < f f \E [f'(B?)f(B?)]\ds 2H dr™ 
Jo Jo Jo Jo 



t rs 



Jo Jo r \ s ~ r ) 
<C H \\f\\ 2 ^ 

for all e > and t > 0. Thus, we get 



\E [f{Bf)f{B?)\] ds 2H dr 2H 



1 

AH 



<C H \\f\\le- 



[ [ ^ £ (s,r,l)ds 2H dr 2 
Jo Jo 

Similarly, we can also obtain the estimate (j4.7j) for k = 2, and the lemma follows. 
Recently, Gradinaru-Nourdin [12] introduced the following perfect result: 



(4.11) 



n s+e 



Y s f(- 



B 



H 



-)ds — >■ E[f(N)} / Y s ds 



□ 



Theorem A (Theorem 2.1 in Gradinaru-Nourdin [12J). Assume that H E (0,1). Let 
f : M — > R be a function satisfying 

(4.10) \f(x)-f(y)\<C\x-y\ a (l + x 2 + y 2 ) b , (C > 0, < a < 1, b > 0), 

for all x,y 6K, and let {Yj : t > 0} be a continuous stochastic process. Then, as e — > 0, 



almost surely, uniformly in t on each compact interval, where N is a standard Gaussian 
random variable. 



According to the theorem above we get the next lemma. 
Lemma 4.2. Let < H < 1 and f £ C(R). We then have 

(4-12) lim * f f{B?)(B? +e - B«) 2 ds 2H = I 

elO £* n J J 

almost surely, for all t > 0. 

As a direct consequence of Lemma 14.21 for / £ C 1 (R) we have 

(4.13) } = 2H f f(B?)s 2H -'ds 

Jo 



H\j2H 



f{B?)ds 
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for all < H < 1. In fact, the Holder continuity of fBm B H yields 

lim * f o(B» +£ - B? )(Bf +e - Bf fds 2H = 
eiO e Jo 

almost surely. It follows that 

= lim ^ /** f(B? )« e - Bf ) V* = /* /(Sf )cfa M 



almost surely. 

Now we can show our main result. 



Proof of Theorem \4 ■ 1\ Given / £ If / € C 1 (M), then the theorem follows from the 
identity (|4.13|) and the follows estimate: 



(srf^dsdr 



<Ch/ / — ^ [/W )] 

jo jo r 2 fs — r 2fl 



<C„ f s^ElfiBf)] ds<C H 
Jo 



2 



by Lemma [331 and (gTUJ) . Let now / Cj 
Consider the function ( on R by 



■5c / 



(4.14) ((x) :-- 



ce(— D 2 - 1 , xG(0,2), 
0, otherwise, 



where c is a normalizing constant such that J K ((x)dx = 1. Define the so-called mollifiers 

(4.15) ( n (x) := n((nx), n = l,2,... 
and the sequence of smooth functions 

(4.16) f n (x)= [ f(x-y)Uy)dy= I f(x-^-)ay)dy, n = l,2,... 

Jr Jo n 

for all x G R. Then {/„} C C°°(M) n Jif and f n converges to / in Jf, as n tends to 
infinity. 

On the other hand, by Lemma 14.11 we have 



P(\Je 1 (f,t)-J £2 (f,t)\>6)<P[\J £l (f-f n ,t)\>-) +P[\J £2 (f-f n ,t)\ > - 

5 



+ Pi |J £1 (f n ,t)- J 6a (/ n) t)| > 3 
< ^ll/-/n|llr + ^(l^i(/n,i)-4 2 (/n,i)| > 
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for all n and 5, £1,62 > 0. Combining this with 



lim J e (J n ,t) = [f n {B H ),B H }f ] =2H [ f n (B^)s 2H - l ds, n > 1 

in probability, we show that the generalized quadratic covariation [f(B H ), B H ]( W ^ exists 
for / G Jtf. Thus, the estimate (j4.4p follows from Lemma 14.11 This completes the 
proof. □ 

Corollary 4.1. Let /, /1, /2, • • • G Jtif- If f n f in , as n tends to infinity, then we 
have 

[f n (B H ),B H ]^ ^[fiB^B 1 

in L 2 as n — )• 00. 

Proof. The convergence follows from 



E 



[f n {B H ),B H ] { P - [f{B H ),B H ] { P 2 < C H \\f n - f\\%, -> 0, 



as n tends to infinity. □ 

By using the above result, we immediately get an extension of ltd formula stated as 
follows. 

Theorem 4.2. Let < H < ^ and let f G Jif be left continuous. If F is an absolutely 
continuous function with the derivative F' = f, then the following ltd type formula holds: 



(4.17) F(B H ) = F(0) + f f{B?)dB* + ~ [f(B H ),B 



t 



Clearly, this is an analogue of Follmer-Protter-Shiryayev's formula (see Eisenbaum [7J, 
Follmer et al [TO], Moret-Nualart [20], Russo-Vallois [26], and the references therein). It 
is an improvement in terms of the hypothesis on / and it is also quite interesting itself. 

Proof of Theorem \4-S?\ If F G C 2 (IR), then this is Ito's formula since 



[f(B*),B H }W=2H [ f'(B?)s 2 ^ds. 

Jo 



For F $l C 2 (M) , by a localization argument we may assume that the function / is uniformly 
bounded. In fact, for any k > we may consider the set 

n k = J sup \Bf I < k \ 

[0<t<T J 

and let f^ be a measurable function such that /t fc l = f on [— k,k] and such that 
vanishes outside. Then is uniformly bounded and f^ G Jif for every k > 0. Set 
jL-F^ = and = F on [—k,k]. If the theorem is true for all uniformly bounded 
functions on Jtf, then we get the desired formula 

(W) 



FW(B? ) = FM(0) + / f [k] {Bf)dB* + - \fW(B H ), B H 
Jo 2 L 



on the set fifc. Letting k tend to infinity we deduce the ltd formula (|4.17|) for all / G Jf 
being left continuous and locally bounded. 
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Let now F' = f E be uniformly bounded and left continuous. For any positive 
integer n we define 

F n (x) := [ F(x - y)Cn(y)dy, 

Jr 

where £ n , n > 1 are the mollifiers defined by (|4.15p . Then F n £ C°°(1R) for all n > 1 and 
the Ito formula 

(4.18) F n (B?) = F n (0) + f f n (B?)dB? + H f f n (B^)s 2H - l ds 

Jo Jo 

holds for all n > 1, where / n = F' n . Moreover using Lebesgue's dominated convergence 

theorem, one can prove that as n — > oo, for each x, 

F n (x) — > F(x), /*(*)— >/(x), 

and {/ n } C f n —tf in as n tends to infinity. It follows that 

jo 

in L 2 by Corollary 14.11 as n tends to infinity. It follows that 

J% n (Bf)dB? = F n (B t H )-F n (0) - ±[f n (B H ),B H ][ W) 

^F{B?)-F(Q)- l ~[f{B H ),B H ]W 

in L 2 , as n tends to infinity. This completes the proof since the integral is closed in I? . □ 

5. Integration with respect to the local time 
In this section we assume that < H < \ and study the integral 

(5.1) / f(x)J? H (dx,t), 

Jr 

where / is a deterministic function and 



Sf H (x,t) = 2H [ 8{Bf -x) 
Jo 



s^ds 



is the weighted local time of fBm B H . Recall that the quadratic covariation [f(B),B] of 
Brownian motion B can be characterized as 



[f(B),B] t = - [ f(x)J? B (dx,t), 
Jr. 



where / is locally square integrable and Jz? s (x,i) is the local time of Brownian motion. 
This is called the Bouleau-Yor identity. More works for this can be found in Bouleau- 
Yor g], Eisenbaum [7J, Follmer et al [ID], Feng-Zhao [9], Peskir [23], Rogers-Walsh [25] . 
Yang-Yan [29], and the references therein. However, this is not true for fractional Brow- 
nian motion. For i < H < 1, Yan et al \2ti\ [27] obtained the following Bouleau-Yor 
identity: 

[f(B H ),B H ]i W) = - [ f(x)J? H (dx,t). 
Jr 

In this section we show that the identity above also holds for < H < i. 
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Take F{x) = (x — a) + — (x — b) + . Then F is absolutely continuous with the derivative 
F' = l(o,6] ^ being left continuous and bounded, and the ltd formula ()4.17j) yields 

[l {aM {B H ),B H ] { p = 2F(B?) - 2F(0) - 2 f \ {aM {Bf)dBf 

J 

= ^ H {a,t)-^ H (b,t) 

for all t G [0, 1]. Thus, the linearity property of generalized quadratic covariation deduces 
the following result. 

Lemma 5.1. For any f A (x) = J2j fA^-ua^x) € &> the integral 

I fA(x)J? H (dx,t) := V/, [^(aj^-^iaj^t)] 
Jr i 



t 



exists and 

(5.2) [ f A (x)J? H (dx,t) = -[f A (B H ),B 

Jr 

for all t G [0, 1]. 

Thanks to the density of $ in 3^ , we can then extend the definition of integration with 
respect to x i— > (x,t) to the elements of J4f in the following manner: 

/ f{x)J? H {dx,t):= lim / f A>n {x)J? H (dx,t) 
Jr n - > °° Jr 

in L 2 for / G J2f provided /A, n - ^ / i n \ as n tends to infinity, where {/A,n} C The 
limit obtained does not depend on the choice of the sequence {/A,n} and represents the 
integral of / with respect to «£f . Together this and Cor pilar v 14 . 1 1 lead to the Bouleau-Yor 
identity 

(5.3) [f(B H ),B H ] ( p = - [ f(x)J? H (dx,t) 
for all t G [0,1]. 

Corollary 5.1. Let < H < ~ and Zei /, /i, /2, . . . G ^\ If fn~* f i n as n tends to 
infinity, we then have 

[ f n (x)J? H (dx,t)^ [ f(x)J? H (dx,t) 

Jr Jr 
in I? , as n tends to infinity. 

According to Theorem 14.21 we get an analogue of Bouleau-Yor 's formula. 

Corollary 5.2. Let < H < ^ and let f G Jif be left continuous. If F is an absolutely 
continuous function with the derivative F' = f , then the following ltd type formula holds: 

(5.4) F(B? ) = F(0) + f f(B?)dB? - \ [ f(x)3? H (dx,t). 

Jo z Jr 



GENERALIZED QUADRATIC COVARIATION 19 

Recall that if F is the difference of two convex functions, then F is an absolutely 
continuous function with derivative of bounded variation. Thus, the Ito-Tanaka formula 

F(B t H ) = F(0)+ f F'(B*)dB? + l I J? H (x,t)F"(dx) 
Jo 1 Jr 

= F(0)+ f F\Bf)dB* - \ [ F'(x)J? H (dx,t) 
Jo 1 Jr 

holds. This is given by Coutin et al [I] (see also Hu et al [17]). 

Remark 1. By the proof similar to Lemma 3.1 in Gradinaru-Nourdin [T2], one can obtain 

the following convergence (see also Gradinaru-Nourdin |13j): 

n ft 
(5.5) lim £ (A,) 2 "- 1 g(B»)(B» - B»_ ± f = / g{B? )s^ds 

j=l J o 

almost surely, where TT n = {0 = to < t\ < ■ ■ ■ < t n = t} denotes an arbitrary partition 
of the interval [0,t] with ||7r n || = supj(tj — tj-i) — > 0, Aj = t - and g G C(M). Thus, 
similar to proof of Theorem 14.11 we can show that the convergence 

n „ 

2H lim {^j) 2H l if(B^) - f(B^ )}(B^ -B»_ 1 ) = - / f(x)J? H (dx,t) 
holds, which deduces 



2H lim Y,( A ^ H ' 1 {f( B ^-f( B lMBS-Bl 1 ) = [f(B H ),B H ]r\ 

3=1 

where / G J(f and the limits are uniform in probability. 



6. The time-dependent case 
In this section we consider the time-dependent case. For a measurable function / on 
R x M + we define the generalized quadratic covariation [f(B H , •), B H y^' of f(B H , •) and 
B H as follows 

(6.1) [f(B H r ),B H )l W) :=hm4r f + e) - /(flf ,*)} « e - Bf )d S 2 * 

for i G [0, T], provided the limit exists uniformly in probability. We prove the existence of 
the quadratic covariation. 

Consider the set Jtff* of measurable functions / on 1 x M + such that the function 
1 1 y f{-,t) is continuous and < +oo, where 



o 



V^s 1 - 11 V Jo Jr ' V2^(T-sy- H 



with ip s (x) = ^ gH e . Then is a Banach space and the set of elementary 
functions of the form 

(6-2) /aOM) =^,fiA(x i - 1 ,x i ](x)l(s j - 1 ,s j }(t) 
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is dense in where {xj,0 < i < n} is an finite sequence of real numbers such that 
Xi < Xj+i, {sj, < j < m} is a subdivision of [0, T] and (/jj) is a matrix of order n x m. 
Moreover, ffl* contains the set 7 with 7 > 2 of measurable functions /onl such that 

|/(x, a )re-5^ ^ dS <oo. 



As a corollary of Theorem A, we have 

(6.3) lim J_ f s™-^g{B»,s)(B^ £ - B»fds = f g(B^s)s 2H - 1 ds 
£ 4-u £ Jo Jo 

almost surely, for all t > if g is continuous. This proves the following identity: 

(6.4) [f(B^,.),B H ]^ =2H J^(Bf,s)s 2H ^ds 

for all t > 0, provided / G C ' (M x R+). Thus, similar to proof of Theorem 14.14 ° ne can 
obtain the next theorem. 

Theorem 6.1. Let < H < \. If f € J^f*, then the generalized quadratic covariation 
[f{B H ,-),B H }^ exists and 



(6.5) E [f(B H ,-),B H ^ {W) 



it 



2 <c H \\f\\% 



for all t G [0,T]. 



By using the above result, we immediately get an extension of Ito formula stated as 
follows. 



Theorem 6.2. Let < H < \ and let F G C 1 ' 1 (M x R + ). Suppose that the function 
J^F = f G Then the ltd type formula 

F(Bf, t) = F(0, 0)+ jf /(Bf , s)dBf + £ ^F(Sf , + ± [/(B* •), B*]™ 
holds. 

Proof. Similar to the proof of Theorem 14.21 we can use smoothing procedure to prove our 
result. The main different key point is the following approximation: 

F n (x,s):= / F(x - y,s - r)( n (y)( n (r)dydr, n ^ 1, 
J Jr 2 

where £ n , n > 1 are the mollifiers defined by (|4.15p , □ 
We next consider the integral 

(6.6) / / f(x,s)^ H {dx,ds), 



where / is a deterministic function. For elementary function f& G S 1 * of the form (|6.2p we 
define integration with respect to local time Jz? H as follows 

f [ fA(x,s)J? H (dx,ds) := J^fij [X H {x u Sj) 
Jo Jr 



'j-l 

(W) 



,/r j (WO 
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for all f G [0,T]. Notice that 

Sf B (x h 8j) ~ 5? H {X U Sj-i) - ^(Xi^Sj) + JSf H (Xi_i, 

= [jSf H (xi,Sj) - JSf H (xi_i, Sj -)] " [^(xi.sj-i) - ^(xi_i, Sj _i)] 
for all i, j. We get the identity 

(6.7) / I f A (x,s)J? H (dx,ds) = - [/ A (B H ,.),B J 
for all i G [0,7"]. Moreover, for / G we can define 

/ / f(x,s)^ H (dx,ds) := lim / [ f An {x,s)^ H {dx,ds), in L 2 
Jo Vr ^°°Jo Vr ' 

for all t G [0, 1] if /A,n — >• / in =^*, as n tends to infinity, where {/A,n} C <^*. 

Theorem 6.3. Zei < -?/ < ^ and / G T/ien the integral (j6.6|) exisis in L 2 and the 
Bouleau- Yor identity takes the form 

(6.8) [f(B H ,.),B H ] ( p =- f f f(x,s)J? H (dx,ds) 

Jo Jr 

for all t G [0,T]. 

Corollary 6.1. Let < H < \, F G C 1 ' 1 (M x R + ) and = f G J%. Then the ltd 
type formula 

F(B H ,t) = F(0,0)+ [ t f(B H ,s)dB H 
Jo 

+ / ^-F(B H ,s)ds-l f [ f(x,s)J? H (dx,ds) 
Jo at 1 Jo Jk 

holds. 

Finally, let us consider the weighted local time of fBm B H with < H < | on a 
continuous curve. Let a(t) denote a continuous function on [0, T]. Then the function 

fa{x,s) = l(-oo,a(s))( x ) 

belongs to J#i, and the integral 

r-t r 

f a (x,s)^ H (dx,ds) 



and the generalized quadratic covariation [f a (B H , •), B H ] ^ W exist in L 2 . By the idea due 
to Eisenbaum [7j and Follmer et al [lQj . as an example, we can show that the process 

f a (x,s)^f H (dx,ds), t>0 
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is increasing and continuous. Thus, we can define the weighted local time of B H with 
< H < \ at a continuous curve t >->■ a(t) by setting 

Sf H (a(-),t)) = f I f a (x,s)J? H (dx,ds). 
Jo Jr 
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